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Abstract. We study the effect of a many-body interaction on inter-band 
oscillations in a two-band Bose-Hubbard model with external Stark force. Weak 
and strong inter-band oscillations are observed, where the latter arise from a 
resonant coupling of the bands. These oscillations collapse and revive due to a 
weak two-body interaction between the atoms. Effective models for oscillations 
in and out of resonance are introduced that provide predictions for the system's 
behaviour, particularly for the time-scales for the collapse and revival of the 
resonant inter-band oscillations. 



Introduction. Recent experiments proved the possibility to study tlie coherent 
dynamics of interacting many-particle systems [l][2l|3]. Such realisations of many-body 
systems with ultra-cold gases in optical lattices have a short but impressive history 
and open immense possibilities for various fields of physics [U |3] . The demonstration 
of the well-known phenomenon of collapse and revival, the latter being a pure quantum 
effect, with ultra-cold atoms bears witness of this coherent evolution of a many-body 
wave function [21 E]. Additionally, the high degree of control in such experiments allows 
a manipulation of many system parameters and makes them particularly interesting 
for various fields of physics as well as future applications [3 \7\. Different ways of 
addressing additional degrees of freedom in such ultra-cold bosonic gases have been 
suggested [5]- 

In the present Fast Track Communication, we discuss a two-band model with an 
additional external force to control the coupling between the two bands. Applying a 
force to atoms in optical lattices leads to Bloch oscillations and is a realization of a 
many-body Wannier-Stark system [S] . The coupling of the low- lying energy bands in 
such systems has been demonstrated in different experiments, e.g. on Landau-Zener 
tunnelling [SJ |H] and the influence of the many-body interaction on a weak coupling of 
the bands has also been studied theoretically [TU]. We go beyond a weak coupling of 
energy bands and consider an isolated two-band system with a strong external force. 
A closed two-band system is an idealisation but can be realised approximatively with 
ultra-cold atoms [J using different techniques as, e.g., super lattices [Ullll]. Besides 
the possibility of experimental realisation, a closed two-band model is also interesting 
as a simple model system. For the latter, we focus on the two lowest energy bands of 
interacting bosons in a optical lattice V{x) — VoCOs{2kLx), with the wave vector of 
the optical lattice = 2'k/Xl. Then, all parameters of the model Hamiltonian just 
depend on this external potential. The parameters can be computed numerically (see, 
e.g.. Appendix A in |13) ) and analytical approximations exist for them for not too 
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small amplitudes Vb Using this setup, we are able to identify regions of strong and 
weak inter-band coupling. A weak two-body interaction introduces new energy scales 
in the coherent evolution of the many-body wave function, leading to a collapse of such 
oscillations, but the form of the interaction gives rise to subsequent revivals. We give 
analytical expressions for all time-scales in this many-body realisation of a two-band 
collapse and revival phenomenon, and experimental ramifications for a realisation with 
ultra-cold bosonic gases are discussed. 

The many-particle model. We study a two-band Bose-Hubbard model with an 
additional external Stark force for a strong coupling of the two bands, introduced in 

[mi] 

n = J2 [^r< - + h.c.) + ^ntint - 1) 

i=\ 

+etn\ + ^±{bl,b, + h.c.) + ^-^n\{n\ - 1) + FC,{b\ai + h.c.) 

+2gW^n'in\ + ?^{b\b\aiai + h.c.) . 

Here, a; {a\) annihilates (creates) a particle at site / of totally L sites in the lower band 
and bi {b\) in the upper band. The corresponding number operators are n° = ajo;, 
n\ = b\bi. The bands are separated by a bandgap A and have onsite-energies 
ej — ±A/2 + IF, respectively. We include hopping between neighbouring sites in 
band a, & with a hopping strength ta,tb > 0, and a repulsive interaction between 
particles occupying the same site in band a (b) with a strength Wa (Wb). The 
two bands are coupled via CqF, with the external Stark force F and a coupling 
constant Co depending on the depth of the lattice Vq [131 E] , and also via the inter- 
particle interaction with a strength Wx- All parameters are measured in recoil energies 
l?roc = fi^fc£/(2m) and we set h — 1 troughout. Focussing on a realisation with a 
single optical lattice rather than a superlattice, the relation between the parameters 
is generally: A 3> ta,tb, as well as ta,tb ~ Wi and Co — 0(10"^). We take the 
external force as a free parameter. We assume that the interaction strength can 
be tuned (experimentally, e.g, by the use of Feshbach resonances [4]) and include 
a corresponding scaling factor g to all interaction terms. To study the occupation 
of the upper band, we prepare the system in an initial state \ip{0)), with a uniform 
distribution of particles in the lower band only and evolve it in time by the many-body 
Schrodinger equation. The quantity we study is the (normalised) number of particles 
in the upper band 

A/;(t)^l(^(t)l5]nf|V(0), (2) 
I 

where N = + ^he total number of particles. We will refer to Nb{t) as 

occupation of the upper band. 

Let us discuss the non- interacting single-particle case "Ho = H(g = 0) first. We 
apply the following transformation involving Bessel functions of the first kind J„(a;) 
to our operators, which is known to remove the hopping terms in the single-band 
case [TS] 

an='^Jl-n{Xa)ai /3„ = ^ J(_„ (xb)^; , (3) 

iez lei 
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with Xi = tijF ^ i — a,b. Using relations for Bessel functions, we arrive at 

no = J2 h""/"; + + CoFJ2 Ji-niAx){al/3i+n + h.c.)] , (4) 

iez n 

where Ax = Xa+xj, and ef = ±A/2 + IF as above. We obtain couphng between any 
two sites of the two different bands, weighted by Bessel functions. In ^ the coupling 
between different and possibly remote sites originates from an on-site coupling and 
subsequent hoppings. It can thus be considered a higher order process in the original 
basis. But in the transformed Hamiltonian (jH), this coupling is now direct with a 
strength modified by the factor J;_„(Ax). 

The Hilbert space in the many-particle problem is spanned by Fock states with 
fixed particle number jn", . . . , n^; nj, . . . , n^) with single-particle basis ai,bi. The total 
dimension of the Fock space for a given number of atoms N and lattice sites L per band 
is given by dimH = {N + 2L — 1)\ / [N\{2L ~ 1)!]. For numerical simulations, we change 
to the interaction-picture with respect to the external force |16j which removes the tilt 
Iri^'^F and replaces a[_|.ia; — e'^*a[_|_j^a; (and likewise for h\j^^bi). The Hamiltonian 
is then time-dependent with a periodicity of Tb = 27: /F and decomposes into a direct 
sum of operators for specific quasi-momenta k [16) . As a consequence, the size of the 
Hilbert space is reduced by a factor of the order of L [T31 [HI US] . Since the different 
subspaces are physically equivalent, we restrain our discussion to the k — 0-subspace 
of the Hilbert space ■ 

For the time-evolution of a given initial state, we use either a direct numerical 
integration (with an adaptive step-size Runge-Kutta algorithm [IT]) or an eigenbasis 
expansion after diagonalising the problem: 

\ip{mTB)) = ^ c„ exp(-ie„mrB)|£„). (5) 

n 

Here we use the eigenstates ?7_F(rB)|e„) — exp(— ie„Ts)|£„) of the Floquet-Bloch 
operator (T denotes time-ordering) 

Uf{Tb) = rexp[ - i / 'H{t)dt\ , (6) 



Jo ^ 

since the Hamiltonian is explicitly time-dependent, with a periodicity Tb- In addition, 
this gives us the full spectrum and enables access to relevant energy scales of the 
problem, as well as an identification of the most important states participating in the 
time-evolution. 

Results. For the specific system under consideration, i.e., bosons in optical 
lattices, both hopping coefficients ta,th are smaller than unity and (since we are 
interested in strong inter-band coupling) we take the external Stark force F to be 
much larger than the hopping coefficients: Ax = {ta+tb)/F <^ 1. The non-interacting 
Hamiltonian Q now allows for simple solutions for two regimes: values of the external 
force F not leading to a degeneracy between energy levels of different bands {off- 
resonant regime) and values of the force F leading to such a degeneracy {resonant 
regime). For the off-resonant regime we make use of Ax ^ 1 and neglect all Bessel 
functions in Q except for Jo(Ax) « 1. The Hamiltonian then decomposes [T5] into 
a sum of independent two-level systems and the occupation of the upper band when 
initially zero follows a simp le Rabi formula Uh{t) = [I + AV(4C^F2)]-i sin^{ht/2), 
where A = -y^ A^ -f ACqF'^. This corresponds to Rabi oscillations between the bands 
with an amplitude much smaller than unity and a period Ta ~ 2t:/A of the order 
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of the Block period Tb = 27: /F. An example is shown in the lower panel of figure [T] 
where the off-resonant contribution to the oscillations is shown. 

Although the coupling from a site I to sites in the other band with different index V 
is usually small (cf. discussion after (HI)), it is important when the two levels become 
degenerate in energy, i.e., for resonant values of the force F. This happens when the 
energy gap between both bands is close to an integer multiple of the external force 
A « rF and we refer to this regime as resonant of order r. In resonance, the coupling 
of the degenerate levels is most important and the Hamiltonian of (HJ can similarly 
be reduced to a sum of independent two-level systems 



and diagonalised by = '^^^^ ^ ^i+r) — '^(l^i " ^i+r)- The resonant 

oscillations (of order r) between the two bands have an amplitude of almost unity 
and a period given by Tj-os — '^/\CoFJr{Ax)\ 3> Tb- An example of these resonant 
oscillations for r = 2 is shown as the thin dashed line in figure [1] The period 
predicted by the reduction to independent two-level systems is T^^"^ = 285 Tb for 
the parameters given there, in excellent agreement with figure [T] where the actual 
period is Tics ~ 288 Tb- We found equally good agreement in numerical simulations 
for different lattice depths and other orders r of resonance not explicitly reported here. 

Let us now study the effect of the many-body interaction in the original 
Hamiltonian ([1]). Figure [T] also shows the occupation of the upper band as a function 
of time for the initial state |?A(0)) = |1, . . . , 1; 0, . . . , 0) in a weakly interacting system 
with g = 0.1. We observe a decay of the resonant oscillations followed by a major 
revival. At later times several minor revivals occur (not shown in figure [1} . This 
effect is stable against variations of the system parameters (as number of particles 
and lattice sites, even for fillings N/L of order but not always close to one), with 
the time scales of decay and revival depending on the interaction strength g (see the 
following section and figure [5]) . We found the same phenomenon with different initial 
states in numerical simulations as long as the particles are mainly delocalised along 
the lattice and occupy only the lower band, i.e., excluding Fock states with all particles 
on one lattice site for instance. 

The oscillatory behaviour depicted in figure [1] strongly reminds of the collapse and 
revival effect known from quantum optics [19_. A specific feature in these systems is 
a linear dependence of the collapse and revival time on the inverse coupling strength, 
i.e. tj-cv oc 1/g and tcoU oc 1/g where g usually denotes the coupling strength between 
the light field and the two- level atom in the quantum optical systems [12] ■ We verified 
numerically for our system that the observed decay and revival times obey a similar 
dependence with the interaction parameter g as coupling strength in our model. We 
define ^coii as the time when the difference between the maximal and average amplitude 
of the inter-band oscillations has fallen to 1 /e, i.e., Afb{tco\\) = 1 /2-I-1/ (2e). The revival 
time is taken as the maximum of the revived oscillations. Figure[2]clearly demonstrates 
the linear dependence of these times on the inverse interaction strength. 

To find an effective description, we try to understand the most relevant interaction 
processes. Clearly, the repulsive two-body interaction disfavours double occupancy of 
sites since this will always cost an interaction energy Wa{b) for two particles occupying 
the same site in the lower (upper) band. Starting from an initial state with population 
of the ground band only, the strong Stark force leads to an occupation of the upper 
band. Doubly occupied sites are also suppressed there, but two particles may sit at 
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Figure 1. (colour online) Occupation of the upper band as a function of time 
for a resonance of order r = 2. Shown are the cases: vanishing two-body 
interaction {g = 0, thin dashed line), weak two-body interaction (g = 0.1) with all 
interaction terms (thick line) and only one interaction term 2Wx Yli '^J'"? (thick 
dashed line) included. The collapse and revival times, tcoii and trcv , are indicated 
by arrows. Lower panel: Magnification of the initial oscillation showing small 
non-resonant oscillations of period on top of the resonant oscillations with 
a much longer period Trcs. Parameters correspond to Vq = 4: A = 4.39, Co = 
-0.15, ta = 0.062, tfc = 0.62,1^0 = 0.030, VFt = 0.018, W^i = 0.012; Af = L = 5 
and F = 2.2207. 




Figure 2. (colour online) Collapse times (open symbols) and revival times (filled 
symbols) versus 1/g for N, L = 5,5 and Vb = 4: □, Vb = 5: O , Vb = 6: 0. Order 
of the resonance; r = 1. The collapse time is defined via A/i,(tcoll) = (1 + e~^)/2 
and the revival time is chosen as the next maximum in the oscillations after initial 
decay t > tcoii, as indicated in figure^ (upper panel). 



the same site in either band, i.e. "on top of each other" . Indeed, the most important 
interaction term in the time-evolution is 2Wx J2i '^""fj since it aheady gives a non- 
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zero contribution when there is only one particle per site in each of the two bands. 
In fact, comparing the time evolution of the initial state with all interaction terms 
and only the one mentioned shows almost no difference (cf. figure [Ij. We focus 
on fillings close to unity h — N/L « 1 and study the time-evolution of states of 
the form j-^o) = |1, 1, . . . , 1; 0, . . . , 0), which is not an eigenstate of the system in 
resonance. Note that this is the most important contribution to the superfluid ground 
state in an expansion in our configuration state Fock basis. For large enough systems 
the superfluid ground state of the untilted system (and the single-occupancy state 
IV'o) likewise) becomes indistinguishable from a coherent state (eq. (66) in [3|) that 
factorises into a product of local coherent states at each site I: 



n(e^"'|vac),) ^^l^;0)'^l¥';0)• 



(8) 



We denote this coherent state with phase = \fn by ly; 0). We are now going to re- 
write this state in the resonant basis and determine the effect of the interaction term 



when acting on this state. We start by inserting the transformation 



li\ = (1/^2) J2n Jl-n{Xa){p^a " ^1) iuto © tO obtain 



V2 



2/2 i 



vac) 



I vac) 



(9) 



= 1 HO] and = 

Let us denote Fock states 
..,n1;n'[,...,nl). From 



where we used J2mez Jm{x)z^"- = exp[a;(z — l/z)/2] for z 
which follows from the properties of the operators a^^\ b^^K 
with single-particle basis ^ii^vi by round parentheses \n^, 
([9]) we see, that the coherent state of the lower band in resonance is a product of 
local coherent states for both bands in the resonance basis Oil'^;"'^)' = 
with (f = ip/\/2. The time-evolution of this state (for the non-interacting system 
in resonance) is simple since it is diagonal in the eigenbasis of the Hamiltonian in 
resonance ([7]) . We continue to study the phase evolution created by the most important 



term 2W^ J^i - 
basis: 



perturbatively by expressing the operators a|t\6m in the fi, u 



exp 



|<p;-<p) = exp 



igW^t 



|^;-(p). (10) 



Since \(p] —ip) is a product of local coherent states we can ignore the sum over / and 
discuss the expected behaviour. Firstly, both Xa,b are much smaller than unity and the 
main contribution in the sums over Zi, . . . , ^4 will come from the zeroth-order Bessel 
functions Jq. Secondly, the product of operators gives 16 different combinations of the 
field operators. But due to the prefactors, the combinations with equal indices are the 
most important. They simply give an integer number when applied to the product of 
local coherent states they are acting on. Taking these two points together, the time 
evolution of this state should show an approximate revival at 

47r 
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This result is valid for large systems and cannot account for the effect of non-universal 
properties like a limited number of particles and lattice sites, but we expect it to yield 
the right order of magnitude for finite systems, and in particular the correct scaling 
with the parameters of external potentials (c.f. figure SI inset, below). 

Additional finite size corrections to (jlip can be understood by using the 
decomposition into the eigenbasis ([S]). The evolution of the occupation of the upper 
band certainly depends on the initial state, which we can take into account by studying 
the weights c„ of the initial state expanded in the eigenbasis. The result of a numerical 
diagonalisation for a system in resonance is depicted in figure [31 where the absolute 
values of the expansion coefficients with their corresponding quasi-energies are shown. 
For vanishing two-body interaction strength g = 0, the quasi-energies from the states 




-0.15 



-0.14 



-0.13 



-0.12 



-0.11 



Figure 3. {colour online) \cn\ versus the corresponding quasi-energies s n 111 fill 
expansion in the eigenbasis of the Floquet-Bloch operator JBjl . Shown are different 
interaction strengths: g = 0.0 (*), g = 0.05 (•), g = 0.1 {■), and g = 0.2 (♦); 
other parameters as in figure [T] We observe that different quasi-energies are 
shifted at different rates when increasing the interaction strength. 



with different occupation numbers are degenerate as expected. The energy difference 
between neighbouring lines of constant quasi-energies corresponds to the time-scale 
of the resonant inter-band oscillations and follows from the diagonalisation of the 
resonant non-interacting system as f^ics — 27r/Tres = 2\CoFJr{Ax)\. In the non- 
interacting system (g = 0), two coefficients are dominating and the difference of the 
quasi-energies yields a single time-scale Tros- When the interaction is turned on, the 
weight of states with many contributions from double- or higher occupancies of sites 
decrease significantly (since they are energetically disfavoured) and their quasi-energies 
are slightly shifted. But, surprisingly, only a limited number of additional coefficients 
Cn contributes significantly in the eigenbasis expansion even for <? ^ 0. The observed 
collapse and revival signal is now determined by a few expansion coefficients that are 
much larger than the others. If we focus on the three largest coefficients, denoted 
as ci,C2,C3 and sorted by their quasi-energies £i,£2,£3, we find that the latter are 
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shifted by the interaction by different amounts. The differences between neighbouring 
quasi-energies, — |e2 — £i \ and W23 — |e3 — £2! (shown for the example g = 0.2 in 
figure [3]), lead to a beating between two oscillations with periods Ti2,T23 « Tics and 
the revival time will thus be given by 

27r T12T23 



(12) 



' W12 123 — T12 

a numerical diagonalisation but gives a clear physical 



W23 

This estimate requires 
interpretation to the revival time observed in a specific realisation with N atoms 
on L lattice sites. Thus, taking (fTT|) and (fT^ together, we have an understanding 
of the general physical mechanism triggering the collapse and revival of the resonant 
two-band oscillations. 

We compare our prediction for the revival time, (fTTT) and (HU, to actual numerical 
simulations in figure 2] We find that pT|) gives the right order of magnitude for 




Figure 4. (colour online) Comparison between estimated revival times according 

to mill ( ) and to II12II (•) with numerical simulations (□) for Vb = 8, r = 1, 

and different system sizes, parametrised by the Hilbert space dimension of the 
K = subspaces. The error bars indicate the width of the revival at half maximum. 
Inset: Scaling of the numerically measured revival time with the lattice depth Vq 
for r = 1. We multiply frev by Wx since Wx also depends on Vq to show the 
remaining non-trivial scaling behaviour. Shown are N,L = 4,4 (dim T-L = 86, 
n), Af,L = 5, 5 (dimW = 402, 0), N,L = &,7 {diraV. = 3876, A) N,L = 7,7 
(dim'H = 11076, x). Filled symbols are frev expected from II12II for the same 
system, and the thick dashed line again shows our universal result l lllll . 



the revival in a specific realisation, and, in particular, it shows the correct scaling 
behaviour with the depth of the optical lattice (shown in the inset of figure [4]). 
Additionally, figure 0] shows that the corrections for specific sizes of the system are well 
accounted for by (jl2p . which only slightly underestimate the revival times by a few 
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percent. This deviation could be corrected by including more than three participating 
states, extending in this manner the arguments which lead to (fT2|) . 
In an experimental realisation, the shorter the collapse time tcoii the easier it could be 
observed, and for an estimate we make use of the fact [121 HI] that the collapse time 
is proportional to the revival time 

^coll — 7T To" ^rcv, (-^^) 

where An denotes the width of the distribution of coefficients c„ necessary to expand 
the initial state in the eigenbasis. For the specific example of Vb = 5, g = 0.1, r = 1 
and N = 5 = L we find An w 0{1) for the width of the distribution (i.e. the g 
couplings include just one or very few additional states as in the derivation of p2p 
above), such that we estimate irov/icoii ~ 3 compared to ircv/icoii ~ 5.7 ±0.1 found 
numerically. 

In this work, we focused on realisations with ultra-cold atoms in a tilted 
periodic potential, and the observed effect can be manipulated by engineering the 
potential [U [TT] or the two-body interaction [3]. Specifically, the revival time dTTI) 
depends sensitively on parameters as the hopping strength and the external force 
close to the zeros of the Besscl function in (jlip . This is analogous to already realised 
manipulations by time-dependent forces as predicted by [22' and realized in [2 3) . The 
observation of Bloch oscillations over thousands of periods and a fine control on the 
two-body interaction have already been demonstrated experimentally Therefore, 
the collapse and revival of resonant inter-band oscillations predicted here should be 
accessible in such state-of-the-art experiments. We finally remark that the collapse 
and revival phenomena discussed in this section have their origin (cf . figure [1]) in the 
degradation (due to interactions) of single particle inter-band oscillations. So, even if 
there are analogies to the collapses and revivals observed in other experiments [2, 6, 24], 
the collapse- revival oscillations reported there arise from the interaction within a 
single-band. Therefore, in contrast to our results, those oscillations would not at 
all occur when the lower band interaction is suppressed, equivalent to Wa = in the 
model here discussed. 

Summary. We studied the coupling between two energy bands in a two-band 
Bose-Hubbard model with an additional tilting force. The force can lead to a strong 
coupling of the bands and we found strong resonances in the inter-band oscillations 
in this lattice model. Furthermore, the two-particle interaction leads to a collapse 
and revival of the resonant inter-band oscillations. Here, we made predictions for the 
relevant time scales which were verified numerically. 

A closed two-band system is an idealisation, but it can approximately be realised in 
various parameter regimes with ultra-cold atoms ^4^ using different techniques as e.g. 
super lattices [TT]. In addition, the use of Feshbach resonances allows a complete 
control and fine tuning of the interaction strength '3| needed to test our predictions. 
Our work is a first important step in going beyond ground-band physics and accessing 
more degrees of freedom in bosonic ultra-cold gases. 
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